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Abstract 

We consider the differential sum rule for the effective scattering rate l/r(w) and optical conduc- 
tivity (Ti{oj) in a dirty BCS superconductor, for arbitrary ratio of the superconducting gap A and 
the normal state constant damping rate l/r. We show that if r is independent of T, the area under 
1/t{oj) does not change between the normal and the superconducting states, i.e., there exists an 
exact differential sum rule for the scattering rate. For any value of the dimensionless parameter 
Ar, the sum rule is exhausted at frequencies controlled by A. We show that in the dirty limit the 
convergence of the differential sum rule for the scattering rate is much faster then the convergence 
of the /—sum rule, but slower then the convergence of the differential sum rule for conductivity. 

PACS numbers: 71.10.Ca,74.20.Fg,74.25.-q 
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Optical constants of solids follow a variety of sum rules. Q| The origins of the sum rules 
can be traced back to fundamental conservation laws and are intimately connected to the 
causality of the electromagnetic response leading to Kramers-Kronig relationships between 
the real and imaginary parts of the optical constants. The analysis of the sum rules is a 
powerful tool to study the distribution of the spectral weight in correlated electron systems. 

The dynamics of conducting carriers is usually described in terms of the effective scat- 
tering rate 1/t{uj) and the effective mass m*(a;) j2| 

^ ^ 47r (T{uoy nil, An u cr{uj)^ ^ ' 

where o{oj) = cri{uj) + ia2{uj) is the complex optical conductivity, = Anne^/mh is the 
plasma frequency, and m^, is the band mass j^]. Under special circumstances (e.g., when 
the Eliashberg theory is valid), the spectra of 1/t{uj) and m* can be linked to the real and 
imaginary parts of the electronic self energy. 

The conductivity, that characterizes the absorbing power of a solid, obeys the /—sum rule 
/o°° dujcTi{(jj) = ujpi/8. This sum rule reflects the conservation of the number of particles and 
physically implies that at a given particle density, the total absorbing power does not depend 
on the details of the interactions, and is determined only by the total number of particles 
in the system It is particularly relevant for a superconductor where the conductivity 
acquires a (5-functional piece proportional to the superfiuid density Ug, and the sum rule 
transforms into J^duai{uj) = UpjS^l — Us/n). As n^/n can be measured independently, 
the sum rule is a valuable tool to study the transformation of the spectral weight associated 
with the superconductivity. 

The issue raised in series of recent works jsl 0, 0, 0] is whether there exists any similar 
sum rule for the relaxation rate. At the first glance, the answer is negative as the integral 
over 1/t{uj) diverges and hence has no physical meaning. This can be readily illustrated with 
a simple Drude model where cr{uj) = {ujpjATi)/{l/T — iuj), i.e., t{ijj) = r, and dLJT^^{Lj) 
is infinite. This implies that there is no sum rule for 1/t{uj) that could be traced to a 
conservation law. However, one still can argue that the frequency integral of the scattering 
rate has some physical meaning. Basov et al. argued, in connection with the cuprates, that 
for experimentally relevant frequencies 1/t{uj) is numerically close to an effective scattering 
rate l/rgr (defined below) that is related to the dielectric function and obeys a physically 
motivated sum rule. In this situation, one can argue that there should exist an approximate 



sum rule if the frequency integral over l/r(a;) is taken in finite limits. This approach was 
thoroughly analyzed by Marsiglio et al. jj]. 

In this paper, we analyze whether there exists an exact sum rule for the full frequency 
integral of the difference between 1/t{uj) at two different temperatures in the normal state, 
and, what is actualy more relevant, between l/tau{uj) in the normal and superconducting 
states. For a Drude metal, the frequency integral over the difference of 1/t{uj) at different 
T in the normal state is either infinite or vanishes depending on whether or not 1/r depends 
on T. We argue that zero and infinity are the only two options for the differential integral 
in any model, and that if t{uj) is T independent, the differential sum rule for the relaxation 
rate is exact between the normal state and a BCS superconductor. 

We also analyze what controls the rate of convergence of the sum rules for conductivity 
and 1/r. This issue is elevant for experimental applications as when sum rules are applied 
to the analysis of the experimental data, one always encounters a problem that actual data 
is available over limited frequency intervals. It is therefore imperative to know rapidly sum 
rules converge to estimate the accuracy of the experiment-based sum rule analysis. 

The analysis in the clean limit Ar ^ 1, where A is the superconducting gap, and r 
is independent on T was performed earlier j^, as a byproduct of the modeling for the 
cuprates. In this limit, A is the largest energy scale in the problem, and the differential 
sum rule between the normal state and a BCS superconductor is exact and is exhausted at 
frequencies of order A. Here we consider arbitrary Ar. We show that even in the dirty limit 
Ar ^ 1, the differential sum rule is exact and is furthermore still exhausted at frequencies 
of order A. This result is not intuitively obvious as the conductivity sum rule in the dirty 
limit is exhausted only at ~ 1/r ^ A. We also found that the functional form of 1/tsc{uj) 
only weakly depends on Ar at uj > 2A. Here and below we use subscripts "sc" and "n" for 
superconducting and normal states, respectively 

The fact that the differential sum rule for 1/r is exausted at oc A in the dirty limit may 
be relevant for the interpretation of the data from some high Tc materials. The frequency 
below which fermionic excitations become incoherent may arise from the energy scale u ~ 
1/r . The convergence of the differential sum rule at 0(A) and the near independence of 
Tsc{uj) on Ar implies that as long as the system possesses a sharp superconducting gap the 
differential sum rule for 1/r is not affected by the increased fermionic incoherence. Indeed 
in high Tc superconductors, 1/r is temperature dependent. The issue of whether or not the 
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differential sum rule is exact between the normal and superconducting states could then 
only be experimentally addressed if measurements could be done in both states at the same 
T. This is obviously impossible. What can be clearly observed, at least in YBCO, is the 
re-distribution of the spectral weight between the normal and superconducting states. If this 
redistribution is over scales of 0(A) then the frequency integral can be restricted to a few 
times A. The approximate sum rule can then be analyzed in the hope that the T dependence 
of the relaxation rate is a minor effect compared to the huge change in 1t{uj) which takes 
place between the normal and superconducting states Note in this regard that since 1/r is 
expressed via both ^(t{uj) and Qa{uj), the analysis of the differential sum rule for 1/r, even 
though it is only approximate, still yields information about the spectral weight distribution 
in a superconductor, which is complimentary to the information from the /—sum rule that 
involves only ^a{uj) 

The issue whether or not 1/r differential sum rule is satisfied in dirty BCS superconduc- 

n 

tors was a subject of recent controversy. Basov et al. p|| argued that in both clean and dirty 
limits, 1/tsc{uj) vanish below 2A, but overshoots the normal state 1/r at larger frequencies. 
In that paper, the profile of 1/tsc{uj) was visually related to that of the fermionic density 
of states, for which the sum rule is exact and reflects the conservation of the number of 
particles. Homes et al. i6| correctly observed that the relaxation rate and conductivity are 
expressed via a current-current correlator, and therefore scale with the joint density of states 
of two fermions about which, they conjectured, no rigorous statements can be made. They 
computed the frequency integral of l/rsc(co') — l/r„(c<j) numerically for a BCS superconductor 
Tn{u!) = T with rA = 1/2 and found that even at u as high as 15A, it is still about 25% 
of its maximum value at = 2 A. They concluded, based on this numerics, that the value 
of sum rule integral depends on the choice for cutoff frequency, and that if the integral is 
truncated at u = 0(A), a finite value is expected. Below we demonstrate explicitly that 
the differential sum rule is in fact exhausted at the energy scale controlled by A, however, a 
rather slow convergence found in Q is real and is due to a weak {\oguj)/u'^ decay of 1/tsc{uj) 
at large frequencies. 

We begin with the general argument that zero and infinity are the only two options for 
doj /t{ijj). This conjecture can be verified by applying the Kubo formula that relates 
cr(co') with the full retarded current-current correlator n(u;): cr(co') = (u;p^/47r) / {—ioj + 
5). Substituting this relation into 1/t{uj) we find that l/r(co') = ^S{uj)/uj where S{uj) = 
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— c<j^/n(c<j). Both ^Jn(c(j) and S{uj) are odd functions of frequency. Since n(a;) is analytic in 
the upper half-plane of complex u and does not have zeros (that can be checked explicitly, 
see below), S{uj) is also an analytic function in the upper half- plane. The analyticity implies 
that, by Kramers-Kronig relations. 



/o t{uj) Jo uj 2 

where C = if the integral converges, and C = oo if it is not. It is easy to show that 
3fJS'(0) vanishes both in the normal and in the superconducting state. In the normal state, 
^U^uj) oc cu^, 5311(0;) oc c<j at the lowest frequencies, hence S^u) oc u'^ and S{0) = 0. In the 
superconducting state 3fJn(0) is finite, while $511(0) = 0, hence S{uj) again scales as cj^, and 
again S{0) = 0. Hence / duo/Tiuo) = C, i.e., it is either zero or infinite. For non-differential 
sum rule, C = 00 as at large frequencies, n(cij) tends to 1, hence S{uj) ^ cj^, and the integral 
in (j2)) diverges. This confirms that there is no physically motivated sum rule for 1/t{uj). 

There are two ways to improve the situation. First, one can introduce an effective Seff{uj) 
that converges at high frequencies and is close to S{uj) at experimentally relevant frequencies. 
Then one can hope to obtain an experimentally meaningful approximate sum rule. For the 
scattering rate, the natural choice, suggested by Basov et al. 3] is to introduce an effective 
scattering rate 

1 r 1 1 

(3) 
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where eiu) = 1 + 'ima{uj)/uj is the dielectric function. For this scattering rate, Seffiyj] 
1 — uj/{uj + 4:T{ia{uj)) is an analytic function, and at high frequencies it vanishes as Seff{uj 
0{l/uj'^). The Kramers-Kronig relations are then applicable, and 

du rcuj': 



Teffiuj) 



3? 



eiuj 



= —i 



-UJ. 



pi 



(4) 



For a Drude metal, 



1 _ 1 ^Pl (r) 

At frequencies u ^ Upi, the correction term in the denominator can be neglected, and 
1/T(.ff{uj) ^ 1/r. As the plasma frequency can well be larger than the fermionic bandwidth, 
the actual integration of experimentally measured 1/t{uj) may not extend to ~ Upi, i.e., in 
the measured range 1/T(.ff{uj) and 1/t{uj) are nearly the same. Still, however, this does not 
imply that there is an "approximate" sum rule for 1/r as / duj/T{uj) and / duj/Teff{uj) both 
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diverge when the integration is restricted to ^ ujpi. Only when the frequency integration 
extends to u > Upi, Jq dx/Teff{x) converges to the sum rule value. In other words, typical 
frequencies for 1/Teff{uj) are of order ujpi, and at these frequencies 1/r and l/T^ff are 
very different. 

As we pointed out above, in this paper we use another approach and analyze the difference 
between 1/t{uj) in the normal and superconducting states. The leading divergent term in 
S{uj) ~ uj"^ at high frequencies is the same for both states. It therefore cancels out in 
Ssc{^) — Sn{uj). If the subleading terms scale as negative powers of frequency at large u, 
Ssc{^^) — Sn{uj) converges and then the sum rule becomes exact. We emphasize that contrary 
to Seffiuj), the difference of Ssd^^) — Sn{uj), if converges, begins falling off at frequencies that 
are still much smaller than the fermionic bandwidth and Upi. In other words, the differential 
sum rule for the scattering rate is exact in a continuous model. 

We now check whether the differential sum rule for 1/r is satisfied in a BCS supercon- 
ductor (we argue that it is exact) and also examine the energy scale at which this sum rule 
is exhausted. In order to explore the issue of convergence we will analyze the differential 
frequency sums for both l/r^u) and ai{uj) defined as 



where Tn = t and, as we pointed out above, the notations "sc" and "n" refer to supercon- 
ducting and normal states, respectively. 

The expression for the current-current polarization operator in a dirty BCS superconduc- 
tor has the form 



where Q± = Q ± uj/2. In the normal state, this reduces to a conventional Drude form 
n„(ti;) = uj/{lj + i/r). In the superconducting state, one can show quite generally that 
53n(u;) vanishes below 2A. At high frequencies, n(co') gradually approaches n(cx3) = 1. 

We first verify whether Sgdyj) — Sniyj) vanishes at = oo. As typical internal frequencies 
in Eq. (jH} are of the same order as external a;, at cx; ^ A, one can expand the integrand in 




(6) 



(7) 
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(jS} in A/cu. The expansion can be straightforwardly carried out for arbitrary At, and the 
result is that at high frequencies 



m.M=mjLu) 1 



1 2AMogf \ 

^ l + (cjr)2 ) 



(9) 



^ (a;r)2 + l lj^ J 

where 3f?n„(cj) = (cjr)^/(l + (ujtY), 5>n„(cj) = — co'r/(l + (tur)^). Substituting these results 
into S{uj) = uj'^/Il{uj) we find that for large u ^ 1/t 



^[Ssc{uj) - Sn{uj)] OC iO ^ '^[Ssc{uj) - Sn{uj)] OC UJ . 



(10) 



We see that both real and imaginary parts of Ssd^^) — Sn{uj) vanish at infinite frequency. 
This implies that Kramers- Kronig transformation is applicable, and hence the differential 
sum rule is exact for a dirty BCS superconductor. 

We now address the issue of the energy scale over which the sum rule is exhausted. We 
consider clean and dirty limits separately. In the clean limit, the frequency integral in 
was evaluated in Ref. 



|. To first order in 1/Ar we have for oj > 2 A 

1 „ / r 4A^^ 



-E 



1- 



TUJ 



(11) 



where E{x) = /[ 



7r/2 



sin^ is the complete Elliptic integral of the second kind 



(note that definitions of E differ in different handbooks). In the two limits E{0) = tt/2 and 
E{1) = 1. The result for ^Ugduj) ~ 1 is valid outside a tiny 0(1/Ar) range near 2A where 
^Usci^^) diverges logarithmically. Substituting Usduj) from (fTTj) into 1/tsc{uj) we obtain for 
> 2A 



1 



1 



E 



•1- 



4A2 



ai{uj) 



1 



(12) 



and l/r(ci;) = cr(to') = for uj < 2A. We plot these functions in Fig. ^ For the differential 
sum rule we then obtain (a; > 2A) 



Uoo) 



2A 



r Jo 



dx 



2A /•! . ^ 
+ / dxE (Vl 



47r2Ar 

We explicitly verified that /r(oo) = 0, i.e., the differential sum rule is indeed satisfied. 



(13) 
(14) 




FIG. 1: (a) The frequency dependence of the relaxation rate t/tsc{'^) in the clean limit, (b) The 
behavior of Irii^) and Ia{^) normalized to their maximum values. 



We also see from ()13|) that at finite tu, It{'^) depends only on cj/A. This implies that in 
the clean limit, the differential sum rule is exhausted at frequencies 0(A). An important 
issue, however, is how rapidly liu) converges to zero at u; ^ A. In Fig. [TJd we plot Iriyj) 
and /cr(i^) evaluated numerically from (jl3j) and normalized to their values at the maximum 
at = 2A. We see that /r(co') converges much more slowly that lai^^)- In particular, at 
uj = 15 A, |/t(^^)| is still about 25% of its maximum value in accord W.K 0. On ..e 
contrary iJuj) is vanishingly small ai uj = 15A. 

valne. T J result fuH, agrees wit. ^. For tire saure pr.t.ca,, vauisK.. 
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A weaker convergence of It{uj) can be understood analytically. Indeed, at high fre- 
quencies, the elliptic function can be expanded in l/u^. This yields 1/tsc{uj) = (l/r)(l + 
2(A/l<j)^ log(2u;/A)). Integrating this expression over frequency, we obtain that at a; ^ A 

/.M = — -[l/2 + log(2.;/A)] (15) 

r LJ 

The conductivity integral meanwhile converges to zero as 

Comparing (|T5|) and ()16|1 . we see that Icr{uj) has an extra (2A/co')^ that accounts for much 
faster convergence of Ia{uj) than of It{^^)- 

We next proceed to the dirty limit Ar ^1. As we said in the introduction, the key issue 
here is whether the sum rule is still exhausted at frequencies 0(A), or one needs to extend 
the integration to frequencies of order 1/r, where the /—sum rule for the conductivity is 
exhausted. At the first glance, in the dirty limit, the frequency integration has to be extended 
to larger frequencies than in the clean limit, as one can easily show that at Ar ^ 1, the jump 
of 1/tsc{uj) at 2 A is small, of order Ar, and therefore 1/tsc{(^) in a superconductor does not 
overshoot the normal state 1/r immediately above 2A. However, as we now demonstrate, 
typical frequencies for the differential sum rule still scale with A. 

The first indication that the physics is still confined to frequencies 0(A) comes from the 
analysis of the form of 1/tsc{uj) at a; ~ r^^ ^ A. Using Eqs. © for n(u;) and substituting 
them into 1/tsc{u)), we obtain 

We see that at frequencies comparable to 1/r, l/rsc{uj) exceeds 1/r, i.e., the overshoot 
occurs at a lower frequency. Further, integrating 1/tsc{uj) from 0(l/r) to infinity we find 
that the contribution to Ir{uj) from these frequencies is of order A^| log(Ar)|. Meanwhile, 
the loss of Ir in the superconducting state between and 2 A is 2A/r that in the dirty limit 
is much larger that A^|log(Ar)|. This implies that even in the dirty limit, the dominant 
contribution to the sum rule comes from frequencies well below 1/r. At these frequencies, 
the current-current polarization operator can be evaluated exactly to leading order in Ar 
as one can pull 1/r out of denominator of the integral in the r.h.s. of (jH)). The remaining 
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FIG. 2: Same as in Fig ^ but in the dirty limit At ^ 1. 



integral is evaluated easily, and on substituting the result into 1/tsc{u!) we obtain for u > 2 A 



1 



UJ 



Tsciuj) r Vcu2 - 4A2 

,2 I 



3? 



UJ 



ai[uj) 



r3? 



E 



4A2 



(18) 



47r V " 4A2 

and l/rsc(ct;) = (Ti(c<j) = for uj < 2A. We plot these functions in Fig|21 The behavior 
of 1/ts(.{oj) near 2A and at high frequencies can be well understood analytically. Near 2A, 
expanding E{x) for large value of the argument, we immediately obtain that 1/tsc{oj) evolves 
continuously (up to corrections 0(Ar)), and very near 2A behaves as 1/tsc{oj) = {1/t)[uj'^ — 
4A2]/ti;2. Still, it overshoots 1/r at u ~ 2.68A and develops a maximum at = 3. 48 A. 
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At larger frequencies, l/Tsdyj) approaches 1/r as 1/tsc{ijj) = (l/r)(l + 2(A/c<j)^log(2a;/A)), 
i.e., exactly the same way as in the clean limit. Analyzing the frequency integrals in Ir{uj) 
and /o-(^) we immediately make sure that they converge at u comparable to 2 A. 

In Fig. |2l3 we compare the rates of convergence of It{^) and Ia{i^) in the dirty limit. 
We see that, as in the clean limit, lai^^) converges better. This time, it is related to the 
presence of the extra logarithmical term in the high frequency expansion of It{uj)- We also 
see that the rate of convergence is almost the same in both clean and dirty limits. We recall 
that this result is not obvious as in the dirty limit, 1/tsc{uj) and cri^u) gradually increases 
above at 2A, while in the clean limit, they both jump at 2A and immediately overshoot the 
normal state values of 1/r and di, respectively. 

To conclude, in this paper we analyzed the differential sum rule for the scattering rate 
and optical conductivity in a dirty BCS superconductor. We demonstrated that this sum 
rule is exact if the normal state r is independent on temperature. For arbitrary At, the sum 
rule is exhausted at frequencies controlled by A, but the convergence is rather weak due to 
logarithmical terms. We showed that in the dirty limit the convergence of the differential 
sum rule for the scattering rate is much faster then the convergence of the /—sum rule for 
the conductivity, but slower than for the differential sum rule for conductivity. The latter 
has the fastest convergence in both clean and dirty limits. 
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